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ABSTRACT 



We study kinetics of single species reactions {"A + A —>■ 0") for general local reactivity Q 
and dynamical exponent z (rms displacement Xt ~ t^^^.) For small molecules z — 2, whilst 
z = 4,8 for certain polymer systems. For dimensions d above the critical value dc = z, 
kinetics are always mean field (MF). Below d^, the density Ut initially follows MF decay, 
no — ~ n^Qt. A 2-body diffusion-controlled regime follows for strongly reactive systems 
{Q > Q* ^ tIq''^^^'^) with no — pa nlx^. For Q < Q*, MF kinetics persist, with n^ ~ l/Qt. 
In all cases n^ ~ l/x^ at the longest times. Our analysis avoids decoupling approximations 
by instead postulating weak physically motivated bounds on correlation functions. 



PACS numbers: 

05.40.-|-j (Fluctuation Phenomena, random processes, and Brownian Motion) 

05.70.Ln (Nonequilibrium Thermodynamics, irreversible processes) 

82.35. -|-t (Polymer reactions and Polymerization) 
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The kinetics of reactions between diffusing particles are anomalous in low spatial dimen- 
sions Mean field (MF) theory, according to which the reaction rate is proportional to 
the product of reactant densities, is only valid for dimensions d above a critical value df.. In 
the most fundamental problem of single-species reactions into inert products, A + A — 0, 
ioi d < dc = 2 the density decays asymptotically as l/t'^/^, independently of the magnitude 
of the local chemical reactivity Q and initial reactant density. That is, in low dimensions 
the 1/Qt decay predicted by MF theory is invalid. This is supported by numerical simu- 
lations 0J^,§,|] and renormalization group studies 0. For the case of "infinitely" reactive 
particles for which Q = oo (probability of reaction per collision unity), rigorous bounds 
on the asymptotic decay of density have been derived and in one dimension exact solutions 
exist [piP,p!0|, 1 1 1 . Many other workers have employed approximate methods starting from 



the hierarchy of coupled dynamical equations for correlation functions of all orders which 
is then truncated by expressing higher order correlation functions in terms of lower order 
correlations |lll,|T2|. 

In this work we address the two major aspects of single species reaction kinetics which 
remain poorly understood. (1) Although it is known that at the shortest timescales MF 
kinetics apply ||l3| , p!^ , the crossover from the short time MF kinetics to the asymptotic 



Xjt'^l'^ behavior has not been established. Some groups on the basis of numerical simulations 
[ll5 i p!6 i |T7| have claimed an "intermediate" time regime during which density decays according 
to a non-universal Q-dependent power law. (2) No systematic theory has been able to predict 
reaction kinetics for general Q H Jl3|Jl^ across all time regimes. As a result numerical 
simulations have been fitted to empirical rate laws |jl8|,|l9|. 

The present work will also consider general dynamics. The above results were for reac- 
tions between particles obeying simple Fickian diffusion for which the dynamical exponent 
z describing rms displacement Xt as a function of time, Xt ~ aitjta)^!^ , is equal to 2. Here 
ta is the time corresponding to diffusion distance of order the particle size a. For reactive 
groups attached to long polymer chains z can also be 4 or 8, depending on time and degree 
of entanglement |^|. In these cases the critical dimension vs, dc = z ||21] below which the 
asymptotic density decay is p2| . 

In this letter we give a complete description of the sequence of kinetic regimes for the 
single species reaction kinetics for arbitrary values of dynamical exponent z (where z is 
assumed a fixed number independent of dimension). Results are derived after postulating 
simple physical bounds on correlation functions without the need to resort to ad-hoc decou- 
pling approximations. In agreement with previous studies, we find that below the critical 
dimension dc = z the density nt decays as at long times. We find two possible kinetic 
sequences for d < dc, depending on the magnitude of Q with respect to a marginal value Q*. 
For Q < Q* ("Weak" systems) the crossover to ~ occurs after a l/t MF regime. For 
Q > Q* { "Strong" systems) the crossover to 1/xf occurs before the 1/t regime has developed. 
In this case a short time regime exists during which the number of reactions is proportional 
to the reactant exploration volume, no — ~ n^xf. 

We begin by noting that the average reaction rate rit = -^Ut at point ri is proportional 
231 to to the 2-body correlation function pt(ri, r2) evaluated at ri = r2. Due to translational 
invariance nt = —Xpt{0, 0) where A = Qa'^. Now the dynamical equation for p involves the 
3-body correlation function p^^\ Following Doi P3[ the exact equation is 
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Pt(ri, Y2)=nl-\ j dr[ dr'^ dt' Gt-t>{ri, ra, r'^, r'^) ptir[, r2)6{r[ - r'^) 

- a| dr[dT',dr',J^dt'Gt-Ar^,r2,r[,r',)pf{T[,r',,r',){6{r[~r',) + S{r'^^ , 

(1) 

where Gt{ri, r2, r'^^, is the net weighting for two particles to arrive at ri, r2 given starting 
points r[,r2, in the absence of reactions. The sink terms on the right hand side of eq. (1) 
describe the three ways in which reactions diminish pt from its initial value Uq. The first 
two-body sink term subtracts off pairs which failed to reach ri,r2 because their members 
reacted with one another at r'^^ at time t'. The remaining two sink terms subtract off pairs 
which would be at ri, r2 but only one member of which reacted at time t' at location rg. Such 
a reaction involves a third particle, weighted by the appropriate 3-body correlation function. 
In eq. (1) we used a 6 function as a reactive sink which is a coarse-grained description of 
the reaction process over a scale of order the particle size a. 

Notice that eq. (1) is not closed in terms of pt since it involves the unknown pi . It is 
in fact impossible to write a closed exact equation for pt since correlation functions of all 



orders are coupled in an infinite hierarchy of dynamical equations [Q. This complication 
may be resolved by approximating 3-body correlations as products of lower order correlation 
functions However we are able to write a closed equation for pt after assuming much less 
restrictive bounds on the magnitude of correlation functions. First we transform eq. (1) to 
an expression for the function gf(r,0) = pt{r,0)/n'f. Notice that gf(0,0) = kt/X where the 
rate constant kf is defined by rif = —ktn'f. We find that eq. (1) is transformed into 

qtiv, 0) = 1 - /* dt' GlZ'ir, 0) k, + iftir) , (2) 
Jo 

where Gl'^^{ri, r2) = / dr2 Gt{r[, r[ + ri, r'g, r2 + r2) is the probability density two particles 
are separated by ri at time t given initial separation r2. The function (pt is defined by 

^,iv)^ J^t' ldr'GlZ,{T,T')k,p,iv') , p,{v) = 2{p,(r|0)-pf)(r|0,0)} . (3) 

Here Pf(r|0) = pt{r,0)/nt and pi (r|0,0) = pl (r, 0, 0)/pt(0, 0) are the conditional densities 
at r, given one and two particles, respectively, at the origin. 

Now we would like to solve eq. (0) for qt{0, 0). This requires information on the properties 
of ft{0) which involves unknown 2-body and 3-body conditional densities. To proceed, let us 
now make the following assumption: We assume that the more particles placed at the origin, 
the lower the conditional density. Chemical reactivity can only induce anticorrelations. Thus: 

pP''(r|0,0) < Pt(r|0) < rit (assumption) . (4) 

Eq. implies the following constraints on p: 

< ptir) < 2nt , f'dt'ke > f dt' f dr' ke pt'{r') . (5) 

Jo Jo J 

The first constraint immediately follows from eq. (|]) while the second is obtained by 
requiring that the magnitude of the second term on the right hand side of eq. (§) exceeds 
that of the 3rd term (since qt < 1) and then integrating over r. 



3 



On the strength of the above constraints, eq. (|^), we will argue below that the solution 
of eq. (0) for gt(0,0), but with the term ipt{0) deleted, gives the correct power law solution 
for kt to within a constant prefactor. Expressing qt{0, 0) in terms of kt, setting r = and 
deleting ipt{0), eq. (0) becomes 

kt = \-\ [ dt'St-t'kf , (6) 
Jo 

where we have introduced the return probability St = Gj'^^(0,0) ~ ^/xf. The validity of 
this approximation is justified by self-consistent arguments whose outline is as follows (full 
details will be published elsewhere |^). If one accepts eq. then one obtains a sequence 



of power law regimes in time for kt and rit (see below). Using these solutions in eq. (|^) and in 
the expression for (pt{0) (eq. (|^)), we have determined the function /i™^^(r') which maximizes 
(/9t(0) (for a given time t) subject to the constraints of eq. (^. This in turn implies an upper 
bound ipf^^^{0) on (pt{0). Now for short times, we find that this bound is much less than 
unity: hence </?t(0) may be deleted in eq. (H) without error and even prefactors are expected 
to be correct. Meanwhile, for long times [rit ~ we find ip'^''^^{0) = A, where A is a 

constant of order unity. We then make a 2nd simple assumption: we assume kt remains a 
power law at long times. We have shown that the deletion of (pt{0) from eq. @ is then 
valid insofar as it will generate the correct power laws for kt (albeit with possibly incorrect 
prefactors) . 

It is now straightforward to solve eq. (|^) for kt. For d > z the integral term in eq. (^ 
is negligible and one recovers the expected MF rate constant k ^ X, true for all times. For 
d < z one finds two regimes: for t < ^2 (see definition eq. (^ below) the integral term in 
eq. (H) is small compared to A and can be neglected. For t > the solution is obtained by 
balancing the integral term with A. One finds 



kt 



xf/t {t* « t) < ' ^ 1 a7[4 \n{t/Q] [tl « t) y^-')' ^ ^) 



as can be verified by direct substitution in eq. (H). Substituting kt in rit = —h^t, it is 
easy to show that for low dimensions two possible sequences of reaction kinetics may occur, 
depending on the magnitude of Q: 

{no - Xnlt no - n^xf l/xf {Q > Q*, Strong) 

{d < z) (8) 
no - Xnlt ^ l/Xt ^ l/xf {Q < Q\ Weak) 

where 

la ta Ano 

(9) 

For the marginal case [d = z) one simply replaces xf in eq. (P) by a'^{t/ta) / hi{t/ta). The 
timescales ^2 ^ind ti are then modified to ^2 = tae^/^^^°-\ ti = ta{noa'^)~^ lia.[l/(noa'^)], while 

Q*ta = l/ln{tl/ta). 

These results for d < z have a simple physical interpretation (see fig. 0). Now if one 
were to assume MF theory were valid, i. e. that particles are distributed at all times as in 
equihbrium and hence rit = — An^, the density would decay as no — Xn^t initially, crossing 
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over to 1/At for t ^ t^. Thus would be the timescale for a particle to react with the 
mean reaction field supplied by the other reactants. For dimensions smaller than the critical 
value dc = z, however, reactions induce nonequilibrium correlations and MF kinetics break 
down beyond a certain timescale. Consider a pair which happened to be initially within 
diffusive range, i.e. within xt, as in fig. [^(a). Since each reactant explores a volume of 
order approximately uniformly, the number of collisions AAcoU after t/ta "steps" increases 
roughly as {t/ta){a'^ /xf) ~ t^~'^/^. Thus for d < dc the net reaction probability QtaMcoW 
is an increasing function of time which becomes of order unity at the timescale of sq. 
(H). For t > t2 a depletion hole of size Xt thus develops in the 2-body correlation function, 
invalidating MF theory. 

For d < dc, the sequence of kinetic regimes after ^2 depends on the relative magnitudes 
of ^2 and t*^. Another relevant timescale is ti, namely the time to diffuse a distance of order 
the typical initial particle separation. One can easily show that the magnitude of always 
lies in between those of ^ and ti. The condition t2 = t*^ = ti defines the critical reactivity 
Q* in eq. (^. Case 1: Strongly reactive species, Q > Q* (or equivalently < < ti)- 

<t <ti reactions are due to the few isolated pairs which happened to be initially within 
diffusive range (fig. |l|(a)). Since most of these pairs will have reacted by t, the number of 
reactions per unit volume during this time regime is proportional to the number of pairs per 
unit volume initially within xt-, n^xf. Now for times t > ti, a region of volume xf contained 
initially several particles. Since only of order one of these particles could have survived by 
t, this implies a 1/xf density decay. Notice that tj^ is irrelevant in strong systems. Case 2: 
Weakly reactive species, Q < Q* (or ti < < t^). Now when nonequilibrium correlations 
develop after the 1/At MF kinetics have already developed. Since > th the same 
reasoning as in the strong case for t > ti implies that a direct crossover occurs after t2 to 
1/xf decay. 

To summarize, we argued that the time dependent rate constant kt in single-species reac- 
tion kinetics is given by a Smoluchowski type expression (eq. (^). Equivalent expressions 
have been the starting point of previous works [25,26,22|. Here we have justified eq. (^ 
based on simple bounds on correlation functions. We argued that the solution for the rate 
constant below the critical dimension will then be correct at short times, but there is ev- 
idence that at long times (during the 1/xf density regime) the actual numerical prefactor 
will be larger than that predicted by eq. (|^). Indeed, for infinitely reactive particles in one 
dimension for z = 2, Torney and McConnell have proved that kt as determined by eq. 
(^) is correct at short times, but is smaller by a factor of 7r/2 from its actual value at long 
times. 

Our result for the density decay (eq. (|^)) does not give evidence for a Q-dependent power 
law decay at intermediate times as suggested in refs. p!5|Jl6|Jl7| for z = 2. We suggest these 
workers may have been observing intermediate kinetics between the 1/t and l/t*^/^ regimes 
ofeq. i). 

We remark that the results here have been for general diffusion dynamics described by the 
dynamical exponent z. They apply for both subdiffusion {z > 2) and superdiffusion {z < 2) 
provided that reactants diffuse independently of one another, and that xt and mean collision 
time ta are well-defined. Implicitly, we have assumed that all moments of displacement are 
characterized by the same scale, Xt- Systems in which reactants perform Levy walks, for 
which the rms displacement does not exist |27,^|, are not covered. The results do not 
apply also to systems in which reactants perform random walks for which the probability 
distribution for the time between successive steps has a power law long-time tail and hence 
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ta does not exist p9| , |30| . In order to describe systems of the latter type, eq. (1) should be 
modified to include an extra integration over all possible collision times. 

To specialize to the most widely studied case of small Fickian molecules one sets z = 2. 
For example, at the critical dimension d = 2 we predict a short time diffusion-controlled 
regime for high reactivity {Q > Q*) in which uq — rit ^ n^t/ Int, and we predict that in all 
cases the long time decay is rit ~ Int/t. 

Finally, it is worth discussing how the above results may generalize to the two-species 
reaction problem, A + i? — > 0, for which segregation of reactants into A-rich and B-rich 
regions occurs at long times in dimensions d < 2z. This leads to an asymptotic decay 

~ {rio/xfy^^ MM- -^^^ d < z, assuming that these kinetics apply for times longer than 
a timescale tg, we may estimate the magnitude of tg by assuming that eq. (||) applies for 
t < ts and demanding continuity in rif. Thus one finds ^ ti for strong systems, while for 
weak systems 4 ~ ia[%«°'(<54)^]^^*-''~^^'' which satisfies < ts < t^. This would suggest 
that at least for weak systems the 1/xf regime does not arise in A + i? ^ 0. Indeed, this 
is consistent with the numerical simulations of ref. ^, where no 1/xf regime was found for 
the case z = 2, d = 1. These workers observed that for sufficiently short times the density 
remained close to no. For high Q values a crossover to occurred at t;. For small Q, 
they found that the regime was preceded by a 1/t regime during t*^ < t < t^. (We 
interpret the onset time here to be tj^, not ti/{Qta) as interpreted by the authors of ref. 
^ The numerical values of these two timescales happen to be very close to one another in 
this simulation.) We hope future numerical simulations as well as experiments on reacting 
polymers where novel values of z are realized will test the validity of the above theoretical 
predictions. 

*** 
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Figures 



FIG. 1. Snapshots of particle distributions during various kinetic regimes for o? < z. Rms 
displacement is Xt and I is typical initial separation among reactants. (a) Times t <^ ti. 
Reactions are due to the few isolated pairs which happened to be initially close enough 
such that their exploration volumes (indicated by dashed lines) overlap by t. For strong 
systems most of these pairs will have reacted by leading to tiq — rit ^ ^o^t • (b) The 
situation at time t = ti. A region of size xf contained of order one particle initially. For 
strong systems reaction kinetics are already diffusion-controlled and a crossover occurs 
to nt ~ (c) Weak systems, U <ti t <ti t^. For weak systems MF kinetics persist 

beyond ti. Many reactants now exist within the exploration volume of a given particle. 
Reaction with the mean density field occurs at a timescale t*^ leading to ~ l/Qt for 
tjl^ <^ t <C ^2- (d) For both weak and strong systems, at sufficiently long times only of 
order one particle survives within a volume of size x^, implying rit ~ l/a^j- 
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